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Ground-State Properties of Extended Two-Channel Kondo Model 

Hiroaki Kusunose* 

Department of Physics, Tohoku University, Sendai 980-77 
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Ground-state properties are examined for an extended two-channel Kondo model where the 
Hilbert space of the localized states is extended to include a singlet state in addition to the 
doublet states. By means of zero-th order variational wavefunctions with different symmetries, 
which are associated with the non-Fermi-liquid and the Fermi-liquid ground states, we demon- 
strate that the channel exchange coupling via the localized singlet state stabilizes the Fermi-liquid 
wavefunction. The ground-state phase diagrams, which are in qualitative agreement with the 
previous study performed by Koga and Shiba, are obtained. The comparison to the structure of 
the resultant wavefunctions suggests that a unique non-Fermi-liquid (Fermi-liquid) fixed point 
exists, irrespective of the localized ground state. 

KEYWORDS: quadmpolar Kondo effect, crystalline electric field, extended two-channel Kondo model, non-Fermi 
liquid, Fermi liquid, variational method 
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In a number of heavy- fermion alloys, remarkable devi- 
ations from the Fermi-liquid behavior have been recently 
observed during the course of investigations of exotic 
heavy- fermion compounds. 1-3 ' This non- Fermi- liquid be- 
havior can be initiated by several microscopic origins 
due to, for example, the two-channel quadrupolar Kondo 
effect, 4 -* a distribution of Kondo temperatures 5 - 1 or the 
critical phenomena in the vicinity of a zero-temperature 
quantum phase transition. 6 ~ 9 ' In most of anomalous al- 
loys, intersite correlations or disorder effects are not 
negligible to understand the non-Fermi-liquid behavior; 
however, a pure single-ion mechanism seems to be essen- 
tial in some dilute alloys. 

In particular, U K Thi_ a; Ru2Si2 (a; < 0.07) shows un- 
usual InT dependence in the specific-heat coefficient 7, 
the c-axis magnetic susceptibility \ an d the a-axis resis- 
tivity p. 2 ) A couple of experimental results scaled with 
uranium concentration convinced us that the anomalous 
behaviors are due to an intrinsic single-ion effect, in par- 
ticular to the magnetic two-channel Kondo effect. 

On the other hand, U a; Lai_ 2 ;Ru2Si2 exhibits an or- 
dinary Kondo effect, 10 ' in which an uranium ion sits 
on a crystalline electric field (CEF) similar to that of 
the Th substituted material under tetragonal symme- 
try. The difference between the two systems is the lo- 
cation of the first excited CEF level under the 5/ 2 con- 
figuration. The first singlet excited state is located at 
around 25K above the time-reversal (T t $) ground dou- 
blet in U a; Lai_ a ;Ru2Si2, while it is located at about 10 3 K 
in U a; Th 1 _ ;c Ru2Si 2 . 

In the case of the CEF-singlet ground state, which 
may be realized in the hexagonal compound UPt3, the 
low-energy physics, even for an impurity case, are still 
unclear; in contrast, extensive experimental investiga- 
tions have been performed in order to understand exotic 
behaviors such as the unconventional superconductivity 
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and extremely weak antiferromagnetism. 11 ) 

From a theoretical point of view, several authors have 
investigated the influence of the CEF-singlet state on 
the non-Fermi-liquid fixed point of two-channel Kondo 
model using the numerical renormalization-group (NRG) 
method. 12, 13 ' They revealed that the mixing between the 
two channels via the CEF-singlet state plays an impor- 
tant role for stabilizing the Fermi-liquid fixed point over 
the non- Fermi-liquid one. Although the properties of the 
Fermi-liquid fixed point are mainly examined, it is not 
yet apparent how to stabilize the non-Fermi-liquid fixed 
point. 

In this letter, it is demonstrated that two ground states 
with different symmetries can be obtained on the basis 
of the variational treatment due to a level-crossing be- 
tween them. 14-16 ' This approach is of great advantage 
for understanding how to stabilize the non-Fermi-liquid 
state over the Fermi-liquid state. 

We begin with an extended two-channel Kondo model 
for tetragonal and hexagonal symmetries, in which the 
time- reversal doublet \E±) and the CEF-singlet \B) 
among the J — 4 multiplet are taken into account and 
only the lowest f 1 intermediate state has been used in the 
Schrieffer- Wolff transformation to derive the exchange 
Hamiltonian. 12 ) Within these restrictions, four states, 
which are specified by channel and spin indices, among 
the j — 5/2 partial wave of conduction electrons are cou- 
pled with the internal degrees of freedom at the impurity 
site. 

Here it is convenient to introduce the local operators 
for the spin and the channel of conduction electrons de- 
fined as 



ha k' a' 
km k f m' 



(1) 

(2) 



respectively. Then the extended two-channel Kondo 
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model is given in terms of both operators by 

E + ,B,E_ 



M - 




J (a 1 + a 2 ) 



Kt 



D/2 



where a h annihilates the fc-Bloch state of conduction 
electrons with the spin a (=f, {) in the channel m 
(= 1, 2). The kinetic energy et spreads over — D ~ +D, 
D being the half of a bandwidth. The exchange couplings 
J, K and |A|, which represent the energy splitting be- 
tween the doublet state \E±) and the singlet state \B), 
are measured in units of D. The definition for each state 
in terms of the azimuthal component of total angular 
momentum, J z or j z , is summarized in Table I. 

Table I. Definition for each state to describe the extended two- 
channel Kondo model. \m) represents the azimuthal component 
of total angular momentum. 



State 


Tetragonal 


Hexagonal 


E± 


a|T3> + /3|± 1} 


a|±4> 


+ /3| T 2} 


B 


(|2> + | - 2»/>/2 


(|3> + | 


-3»/V2 


a ki] 


a|5/2}+/3|-3/2> 


|3/2> 


1 - 3/2} 


a kn 


|l/2> 


1 - 1/2) 


|5/2> 


a fc2T 


1 - 1/2) 


|l/2> 


1 ~ 5/2> 


a k2i 


a|-5/2>+/3|3/2> 


1 - 3/2> 


13/2) 



The second term in eq. (3) is reduced to the usual 
two-channel Kondo exchange interaction with the iso- 
lated singlet state \B) in the case of K = 0. Nonzero 
K term plays the role of mixing the channel degrees of 
freedom and the localized doublet state is simultaneously 
transferred to the singlet state and vice versa. The num- 
ber of conduction electrons in each channel is no longer 
conserved, however, the total number of electrons, 

Q = E Wkm^hma ~\\, (4) 
feme- 
is still conserved. 

Since the Hamiltonian (3) is invariant under the trans- 
formation (maE±) <-> (fhaE^), namely, due to the time- 
reversal symmetry in the original representation of total 
angular momentum as given in Table I, we can revive a 
conserved quantity, 12 -* total spin, J z = ^2j z by assigning 
fictitious spin j z to each state as 



(it) 

(2 1) 



(+2). 



(2T)- 
(E±) 



(+5), 
- (±1), 



(U) 

(B)- 



(0). 



(5) 



Now we construct possible variational wavefunctions 
for ground states of the extended two-channel Kondo 
model. A half-filled Q = state with spin polarization 
of conduction electrons around the impurity will be a 
ground state. The non-Fermi-liquid (NFL) ground state 
must include the overscreening state, at at, | + 1) (only 



one of the time-reversal pair is presented), which is in- 
volved in (Q, J z ) = (0, —1) state, \ j z ) being the localized 
state j z and the state filled inside the Fermi sphere to 
the extent kp. While the nondegenerate Fermi- liquid 
(FL) state should be specified by the conserved quanti- 
ties, (0,0). 

The variational wavefunctions can be obtained 14 ) if 
we begin with the states in which two excited electrons 
above the Fermi level couple with the localized spin state 
and then collect higher-order states with particle-hole 
excitation which are connected with the next lower-order 
one by the exchange terms in eq. (1). We can write a 
zeroth order wavefunction |0, — 1) for NFL and |0, 0) for 
FL states in consideration of the time-reversal symmetry, 
J z «— » J z , as 



|o,-i) = 



>fe F 

E 

1,2 l 



a&lfcil T,fc 2 2|;-l)+< 2 ; |fc 1 2 t,fc 2 l 



(2) I 



+ 0$ |fcx2 I, fe 2 l |; +1) + a£> |fei2 |, k 2 2 |; 0) 



(4), 



, (6) 



|0,0) 



>fe F 



E 

1,2 



+ §/3f 2 } ( Ifcil I, fel I; +1) - |fei2 T, k 2 2 T; -1) 



P12 ( l^i 1 1> fe 2i I; -1) - |fei2 1, fe 2 2 T; +1) 



+ /3 1 3 2 ) |fc 1 2 4,fe 2 lt;0) 



(7) 



where we define I fci ma, kim'a';j z ) = at a\ , \j z ) and 
k n is simply expressed as n in the variational coeffi- 
cients, aW and (3^ . The Pauli principle requires that 



-(3^1 and (3^ = P^i is due to time-reversal sym- 



metry. It is noted that |fc x 1 |, fc 2 l |; 0) (|fcil |, fe 2 2 |; 0)) 
does not connect with the above wavefunction 1 0, — 1) 
( 1 0, 0} ) despite the fact that its conserved quantities are 
the same. 

If we ignore the k dependence of the variational wave- 
functions, they are similar to the wavefunctions which 
are obtained by the diagonalization of zeroth-site Hamil- 
tonian in the NRG method. 17 ) The successive renormal- 
ization procedure can improve these localized wavefunc- 
tions to involve the low-energy states of conduction elec- 
trons. In our zeroth variational wavefunctions, such a 
strict improvement is replaced only to take into account 
the k dependence of conduction electrons. 

The candidates in the variational functions are con- 
nected with each other by exchange interactions. The 
relations among them are schematically shown in Fig. 1 . 
The circle with the number i stands for the z-th compo- 
nent of the variational function whose energy is indicated 
above the circle (the kinetic energy ei + e 2 is omitted). 
We can easily expect that J stabilizes the state (0,±1), 
while the state (0, 0) can be stabilized due to the ex- 
change K because the latter state has two mixing paths 
with respect to K, while in the former state, two mixing 
paths are available for J instead of K. Since the relation 
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J 



1 HW 

J J K 




Fig. 1. The relations among each component in the variational 
wavefunction (a) |0, ±1), (b) |0, 0). The circle with the number i 
stands for the i-th component of the variational function whose 
energy is indicated above the circle (the kinetic energy ei + 12 is 
omitted). 



of the number of mixing paths remains the same when 
considering higher order wavefunctions, we expect that 
the zeroth order wavefunctions contain essential features 
to clarify level-crossing of ground states. 

The Schrodinger equation [H - E(Q, J z )) \Q, J z ) = 
provides us with the following set of equations to deter- 
and E: 



mine a 12 



and $> 



J [A. 



t „(2) 7/4(2) .(3) 



(1) 
2-1 

.(2) 



A 



(3) 



C (3) 



J(A 

(4) 



(3) 
2+1 



A 



(i) 



(62 + A) a% - IK A 



1(2) 




(8) 



?12Pl2 
4l2Pi2 



J ( M-l 



12 



J IB, 



(2) 
2+1 



✓5B< 2 > 



- V8KB[ 3) = 0, 



B, 



(i) 



V2K (B™ 



B. 



(i) 



))=0, (9) 



= 0, 



where £12 = ei + £2 — -E 1 , arid for (C, 7) 



*) 



(A, a) or(S,j8) 

3 Im3i 



First, we consider the case J, K ^> D = 1, which 
is related to the local problem in the NRG method as 
mentioned above. We can assume 7^ = 7* (independent 
of k); then eqs. (8) and (9) are reduced to 



/ 


c 


-j 







C J 




—j 


J (-J 


V 





-2K 



c 

-2V2K 



-2V2K 

C + A 





-2K 


C + A / 

P 3 



Q 



V « 4 / 

= 0, 



0, (10) 



(11) 



respectively, where (C + E)/2 = (e k ) = J2k €k / N an( l * ne 
component of f3 2 vanishes because of the antisymmetric 
condition. 

From these equations, we can obtain ground-state en- 
ergies E(Q, J Z )=E - 2{e k ) as 



E(0,±l) = min(-2J, A), 
E(0,0) = min(0, A) 



(12a) 
(12b) 



for J > K and 
E(0,±1) = min 



i (a - Va 2 + 16K 2 ) , — j 



(13a) 



S(0, 0) = i (A - VA 2 + 32if 2 ) (13b) 

for J <C K. Thus, the ground state (0,±1) can be real- 
ized for the case J > K, | A| or A > J, K, while (0, 0) 
for > J, |A| or-A> J,K. 

Next, we investigate the ground state properties for 
various J and K numerically. The variational param- 
eters are adjusted suitably with advantage both in the 
exchange energies J, K and the kinetic energy e^. The 
phase diagrams of the ground state are shown in Fig. 2(a) 
for A = 0.2 (time-reversal-doublet ground state) and 
2(b) for A = -0.2 (CEF-singlet ground state). The 
boundary lines are determined by E(0, 0) = E(0, ±1) 
and the singlet state is stable for £(0,0) < E(0, ±1), 
while the doublet state is stable for E(0,±1) < E(0,0). 
These results are in qualitative agreement with the pre- 
vious study using the NRG calculation. 12 ) The critical 
strength of the exchange coupling K c for (a) or J c for (b) 
is roughly determined by | A | . It is noted that the bound- 
ary for negative A is shifted to the right more than that 
expected as K ~ J for large J, K. Typical examples for 
the state (0, ±1) are given in the regions I and IV and 
for the state (0, 0) in the regions II and III in Fig. 2. 
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/ IV 




/ (0+1) 







0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 0.6 
J J 



Fig. 2. The ground-state diagrams for exchange couplings J and 
K: (a) A = 0.2 (b) A = -0.2. The nature of the typical ground 
states can be seen in regions I— IV. 



The changes of variational parameters at k = k' = fep 
(k — fcp, k' = kp + 5 for f3^ 2 ') across the phase boundary 
are shown in Fig. 3(a) for J = 0.3 and A = 0.2 (I — > II) 
and 3(b) for K = 0.3 and A = -0.2 (III -> IV). 

For typical (0, ±1) ground state in region 
are dominant. The term associated with the over- 
screening state becomes most dominant as expected in 
the conventional two-channel Kondo model. This situa- 
tion is also realized in region IV, although the localized 
singlet state is stable (in this case, term is expected 
to be dominant), because the exchange energy K sur- 
passes the loss of the localized energy |A|. Since the 
distribution of intensity spreads over the higher kinetic- 
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(a) J = 0.3, A = 0.2 (b)K = 0.3, A = -0.2 
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Fig. 3. Changes of variational parameters at k = k' = fcp (k = 
k F , k' = k F + 6 for /?( 2 >) across the phase boundary: (a) for 
J = 0.3 and A = 0.2 and (b) for K = 0.3 and A = -0.2. 



energy states due to the strong coupling constants, each 
intensity of variational parameters at the Fermi level in 
the latter example is smaller than that of the former ex- 
ample. On the other hand, the (3^ and the (3^ terms 
are dominant for (0, 0) ground state in regions II and III, 
which include the decoupled CEF-singlet state as a limit 
^(i) _> 0. 12) This result implies that a unique non- Fermi- 
liquid (Fermi-liquid) fixed point exists, irrespective of the 
CEF ground state. 

The structure of the wavefunction changes drastically 
but continuously across the phase boundary in order to 
make full use of the larger exchange coupling. As a direct 
result of this change, one stable ground state is replaced 
by the other one. The low-lying excitations between both 
sides of the phase boundary might be completely dif- 
ferent because the nature of the ground state changes 
abruptly at the boundary. This situation is essentially 
analogous to that of an anisotropic spin S — 1 in the 
magnetic field, which is described by the Hamiltonian, 
H = DS 2 Z — S Z B. In this example, the ground state |0) 
for B < D is replaced by | + 1) for B > D. However, our 
variational state has a corresponding many-body reso- 
nance in contrast to the purely one-body state \S Z ). 

In summary, we have investigated the ground-state 
properties of the extended two-channel Kondo model on 
the basis of the variational treatment. We construct 
the ground-state wavefunctions with different symme- 
tries, which are associated with the non-Fermi-liquid 
and the Fermi-liquid ground states. Within the 0th or- 
der wavefunctions, in which the particle-hole excitations 
are ignored, the non-Fermi-liquid ground state is real- 
ized due to the spin exchange J, while the channel ex- 



change K stabilizes the Fermi-liquid state both for the 
CEF time-reversal-doublet and the CEF-singlet ground 
states. This quantum phase transition comes from the 
competition between the different types of ground states 
with corresponding many-body resonances. The com- 
parison to the structure of the resultant wavefunctions 
suggests that a unique non-Fermi-liquid (Fermi-liquid) 
fixed point exists, irrespective of the CEF ground state. 
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